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I. INTRODUCTION 


In this contribution we determine the lowest-order Quantum Electrodynamic corresctions 
and other bound state properties in some two-electron ions with relatively large nuclear 
charges Q, where Q < 28. In general, the lowest-order QED corrections are relatively small 
~ for the ground and low-excited states in the two-electron ions. These QED corrections 
arise from interactions of the bound electrons in ions and virtual photons. General theory of 
the lowest-order QED corrections in atoms/ions was considered in |l|. Since the hrst edition 
of that book the problem of lowest-order QED corrections was re-considered a number of 
times in different books and textbooks (see, e.g., j^, j^). Note that often the lowest-order 
QED correction is called the ‘Lamb shift’,or even ‘Bethe logarithm’. In reality, such a 
correction always contains a number of terms and Bethe logarithm is only one of them. Two 
other terms are the electron-electron contact term and Araki-Sucher term which contains the 
expectation value of the singular operator (^). The explicit formula for the lowest-order 

^ee 

QED correction is given in the third Section III. The formula presented there is based on an 
assumption that the nuclear mass is inhnite, but in actual atomic systems the mass of the 
central ‘heavy’ nucleus is always hnite. Therefore, it is important to have a formula which 
allows one to determine the hrst-order correction upon the hnite nuclear mass. 

By using highly accurate wave functions constructed in this study we can evaluate a 
number of other bound state properties for these two-electron ions, including the {r~^) and 
(r“^) expectation values, expectation value of the Vinti quasi-singular operator and 

(^ew) expectation values for fc = 2,4, 6,... which are needed to evaluate angular dependence 
of the elastic scattering cross-sections of the fast electrons by the two-electron, helium-like 
ions in the Born approximation. We also consider an important problem of Q—dependence 
for a number of expectation values determined in this study. 


II. BOUND STATE WAVE FUNCTIONS 


Modern highly accurate calculations of the lowest-order QED corrections in few-electron 
atomic systems, including two-electron ions, are performed with the use of the non- 
relativistic wave functions which are determined as the solutions of the Schrodinger equation 
^ for the bound states = ET, where E < 0 and H is the non-relativistic Hamiltonian 
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of the two-electron ion with the nuclear electric charge Qe: 


H 



+ V2 + 



Qe^ Qe^ 

r32 rsi r2i 


( 1 ) 


where Vj = and i = 1,2, 3(= N), where the notation N (also 3) stands for the 

central atomic nucleus. In this equation and below the notation h stands for the reduced 
Planck constant, i.e. h = ^, and e is the elementary electric charge. Note also that Eq. (fT]) 
and everywhere below in this study the particles 1 and 2 mean the electrons, while the 
particle 3 (or N) is the atomic nucleus with the mass Mjv 3> mg. In real atomic calculations 
such a mass can also be an inhnite. For light atoms and ions it is very convenient to perform 
all bound state calculations in atomic units where h = l,me = 1 and e = 1. In these 
units the velocity of light in vacuum c numerically coincides with the inverse value of the 
dimensionless hne structure constant, i.e. c = a~^, where a = ^ ^ 7.2973525698-10“^ j^. 
In atomic units the Hamiltonian from Eq.([T]) is written in the form 




^31 ^21 


( 2 ) 


where the nuclear mass must be expressed in the electron mass rUe- 

Formally, in calculations of the lowest-order QED corrections one needs to use the truly 
relativistic wave functions determined from the corresponding Dirac equation for the bound 
atomic states. At this moment, however, highly accurate solutions of the Dirac equation 
can be derived for one-electron atoms and ions only. For two-electron atoms and ions the 
relativistic wave functions can be constructed only approximately and overall accuracy of 
such an approximation is not high. It follows from the well known fact that the electron- 
electron repulsion presents a great difficulty for analytical consideration in the Dirac theory 
of two-electron ions/atoms. In actual computations such an accuracy is substantially smaller 
than analogous accuracy obtained with the use of variational wave functions which were de¬ 
termined as highly accurate solutions of the non-relativistic Schrodinger equation. In this 
study to determine some properties of the two-electron atom(s) and ions we also apply the 
non-relativistic wave functions. All wave functions used in our analysis have been obtained 
as the solutions of the non-relativistic Schrodinger equation with the Hamiltonian, Eq.(l2]). 
Note that this approach with the non-relativistic wave functions works, if (and only if) we 
can derive the non-relativistic variational wave functions of very high accuracy. Indeed, all 
expectation values below are determined in respect with the rules and formulas of the per- 


3 







turbation theory. In such calculations it is extremely important to have extremely accurate 
non-relativistic wave functions constructed in zero-order approximation. 

In this study highly accurate wave functions of the ground 1^S'(L = 0)—states of the 
two-electron ions/atoms are constructed in the form of the following variational expansion 
3 , e.g., je] and Q) 

N 

= (l + A 2 ) exp(-air32 - Argi - y^rai) (3) 


see, 


2=1 


T.c. 


exp(-air32 - Pir^i - 7i^2i) + exp(-^ir32 - - 7ir2i) 


2=1 

The expansion, Eq.Q, is called the exponential variational expansion in the relative coordi¬ 
nates r 32 , r 3 i and r 2 i. Each of the three relative coordinates rij is dehned as the difference be¬ 
tween the corresponding Cartesian coordinates of the two particles, e.g., Vij =| ri — rj \= rji. 
It follows from this dehnition that the relative coordinates r 32 ,r 3 i and r 2 i are translation- 
ally and rotationally invariant. The coefficients Ct are the linear variational parameters of 
the expansion, Eq.(l3]), while the parameters ai,/3i and 7 * are the non-linear parameters of 
this expansion. In general, the total energy of the ground S'—state of the two-electron 
ion depends upon the total number of basis functions N, Eq.(|3]), used in calculations. The 
operator P 12 in Eq.Q is the permutation operator for two identical particles (electrons). 
Very high efficiency of the variational expansion, Eq.(j3]), in actual applications to the two- 
electron ions is related to the fact that all non-linear parameters f3i and 7 * are carefully 
varied in such calculations. The overall accuracy of our variational expansion, Eq.(l3]), for 
all two-electron ions mentioned in this study can be considered as very high and sufficient 
for all our current needs. For the neutral helium atom we have performed numerical cal¬ 
culations with an extended numerical precision. This allows us to determine 34 - 35 stable 
decimal digits in the total energy. By performing these calculations for the helium atom we 
investigated numerical stability of our results for the helium atom and other two-electron 
ions mentioned in this study. It was found in these research that our ‘usual’ wave functions 
which provide 23 - 25 stable decimal digits are sufficient for our purposes in this study. 


III. LOWEST-ORDER QED CORRECTIONS 

Numerical calculation of the lowest-order QED correction in few-electron ion/atom is the 
two-stage procedure. At the hrst stage one determines the lowest-order QED correction for 
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the model ion/atom with the inhnitely heavy nucleus, e.g., for the °°Ne^^ and ions. 

Such a correction is designated below as The explicit formula for the lowest-order 

QED correction AE^^^ in the two-electron ion with inhnitely heavy nucleus is written in 
the form (in atomic units) 


= -Qa^ ^-2\na-\nK^ 


,rl9 
L30 

7 3/ 1 V 

evr" ^ 

KJI\ » gg 


(5(r,^)) + a3 ^ In + 1)1 (5(ree)) 

L 15 3 6 J 

( 4 ) 


where a = ^ = 7.2973525698 ■ 10“^ is the hne structure constant (see above), Q is the 
nuclear charge (in atomic units) and S is the total electron spin. The ground states in all 
two-electron ions considered in this study are the singlet states with S' = 0. Also, in this 
formula In Tfo is the Bethe logarithm (see, e.g., Q. B). 

The structure of the three-term formula, Eq.(jl]), has the following physical sense. Interac¬ 
tions of the bound electrons in ions and virtual photons leads to the variations in the electron 
density distribution. To determine the lowest-order QED corrections in two-electron ions the 
crucial factors are the electron density distributions at the atomic nucleus and at the second 
electron. This expalins the appearence of the two contact terms in Eq.(jll): electron-electron 
nucleus contact term (~ (5(reAr)) and electron-electron contact term (~ ((5(ree)). The last 
term in Eq. (jl]) is called the Araki-Sucher term, or Araki-Sucher correction, since this correc¬ 
tion was obtained and investigated for the first time by Araki and Sucher 8|, js]. Formally, 
this term represents the spatial derivative of the electron density at the nearest vicinity of 
the electron-electron collision point. Note that the expectation value of the term (-|-) is 

^ee 

singular, i.e., it contains the regular part and non-zero divergent part which is also called the 
principal part of the expectation value. General theory of the singular exponential integrals 
was developed in our earlier works (see, e.g., and references therein). In particular, in 
|l^ we have shown that the (^) expectation value is determined by the following formula 


(:k> = G)H + M'5(r«)) 


( 8 ) 


where (-1-)^ is the regular part of this expectation value and (5(ree)) is the expectation value 

^ee 

of the electron-electron delta-function. The presence of non-zero divergent (or singular) parts 
in singular expectation values directly follows from the fact that the corresponding operators 
are self-conjugate. It is interesting to note that ‘to be self-conjugate’ is the most fundamental 
and important property of any ‘symmetric’ operator in QED. All other properties of such 
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operators can be considered as secondary. Here we cannot discnss this interesting, bnt non¬ 
trivial problem. Briefly, we can only say that the overall contribntion of the principal (or 
singular) part of the ^ operator is reduced to the expectation value of the electron-electron 
delta-function with an additional ‘spherical’ factor 4 tt. Formally, in [l^ the equality, Eq.(l5]), 
was derived for the exponential variational expansion, Eq.(l3]), only. However, it can be shown 
that the same equality is true in the general case. Analogous formula can be written for the 
electron-nucleus expectation value (-^) which arises below. 

‘’’eN 

For the actual two-electron ions with the finite nuclear masses we need to evaluate the 
corresponding recoil correction to the lowest-order QED correction, Eq.(|l]). In atomic units 
such a correction to the lowest-order QED correction is written in the form 


- ( 


+ 


VM M +1 




4Q2a3 31 2 , 

H--1-In a — 4 In Kq 

3M ^3 Q ■ 


+ 


7a^ 1 


SnM rl^ 


(<^(rejv)) 

( 6 ) 


where M 3> me is the nuclear mass. The difference AE^j^^ — AE^^^ is the recoil shift of 
the lowest order QED correction. In some earlier works the recoil shift was defined as the 
absolute value of this difference and expressed in MegaHertz, or MHz, for short (1 a.u. = 
6.579683920729T0® MHz). The inverse mass in Eq.([6]) is a small parameter which is 


m 


atomic units) smaller than < 2.5 • 10 ° for all isotopes of the two-electron ions which are 
heavier than the nitrogen two-electron ion N®"*". The dimensionless ratio R = 


QED 




AE 


.QED 


is always small (very small) and can be evaluated asi?~-g~l-10 ® — 1-10 for 

two-electron ions considered in this study. 


IV. SPECIFIC COMPONENT OF THE ISOTOPIC SHIFTS. VINTI OPERATOR. 


Numerical calculations of the isotope shifts in few-electron atoms are performed with the 
use of the formula (see, e.g., |l0| ) 


771 _ 771 I Tpnms I TQsms _ f-I \ 77' 1 rpsms 

Em-Eoc + E^ +Ej^ - I i- -—\Eoo + Ej^ Q) 

^ rUe + 

where Em and E^^ are the total energies of the actual ion/atom with the nucleus, which has 
a finite nuclear mass (M), and model ion/atom with the infinitely heavy nucleus. In this 
equation the notations E^^ and Ef^^ stand for the non-specific and specific mass shifts 
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and nie is the electron mass. Physical meaning of the valne is clear from Eq.([7]). The 
explicit formnla for the specihc mass shift in the A-electron atom/ion takes the form 


E 


sms 

M 


h 

MrUe 


(^oo 


A-1 A 


E E v.-v, 

i=l j=2(J>i) 


( 8 ) 


where /i = is the rednced electron mass and M is the nnclear mass, while rrie is 

the electron mass. The index cxo designates the wave fnnction which correspond to the 
A—electron atom/ion with the inhnitely heavy nnclens. In atomic nnits the formula, Eq.ljHi). 
is written in the form 


jpsms 

— 


1 


M + 1 


(^C 


A-1 A 

E E 


Vi ■ Vj I Too) 


(9) 


i=l j=2(J>i) 

Below we replace the notation Too for the wave function by the traditional notation T. This 
simplihes many formulas which follow. 

For two-electron atoms/ions considered in this study the right-hand side of Eq.Q is 
reduced to one expectation value (T | Vr V2 | T) = (T | pi-p2 | T), where p^ = —iVk is the 
electron momentum of the fc—th electron. The expression for this expectation value can be 
expressed in a number of ways, e.g., (Vi-V 2 ) = — 2 |(Vg)-|-|(V^), where —^(Vg) essentially 
coincides with the single-electron kinetic energy, while — |(V^) is the kinetic energy of the 
nucleus. This relations between the (T | pi-p2 | T) expectation value, single-electron kinetic 
energy and kinetic energy of the nucleus follows from the conservation of the momentum of 
the three-particle system, e.g., in the center-of-mass we can write Pat = pi -|- p 2 . It follows 
from this equation that (P^) = 2 (pg) -|- 2 (pi ■ P 2 ). 

However, in the middle of 1930’s it was assumed that computation of the expectation 
value of any differential operator, e.g., pi ■ p 2 and/or Pg cannot give a numerically stable 
result. In that time people tried to express such ‘unstable’ expectation values in terms of 
the linear combinations of other expectation values which are the functions of the electron- 
nucleus and electron-electron variables. In particular, the following formula for the (pi ■ P 2 ) 
expectation value has been derived by Vinti 11| 


/ \ ^/^31 ■ 1 ' 32 ^ 1 / 1 \ 

(Pi-P2) = Q(—^) + -(—) 


( 10 ) 


where Q (or Qe) is the electric charge of the nucleus (particle 3), and particles 1 and 2 are 
the two electrons. Numerical calculations of the expectation values of the electron-electron 
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repulsion {r^i) = {fee) straightforward. However, the hrst expectation value in the right- 
hand side of Eq. llTnll contains a singular expression. Indeed, in the relative coordinates the 
expectation value ((rai ■ r 32 )/r|^) is represented in the form 


,r3i ■ r32 




' 2V3i'^2^V?/ '■ 


21 


( 11 ) 


31 ^ ' 31 ^ ' 31 '31 

where the last two expressions in the right-hand side are singular, if they are considered sepa¬ 
rately of each other. For years this fact made impossible any highly accurate computation of 
the expectation values of the Vinti operator • Recent analysis, however, indicates that 

’’31 

computations of this expectation value in the exponential expansion, Eq. flTTl) . is reduced to 
operations with the Frullanian integrals only. By the Frullanian (three-particle) integrals 
we mean three-particle integrals which look like singular, but, in reality, they are regular 
integrals 12|], 13|. The expectation value of the Vinti operator contains singularities in 
the second and third terms in the right-hand side, but these singularities cancell each other 
from the hnal expression which becomes regular (for more details, see, e.g., ji^h Note that 
in atcual computations it is better to consider the expectation value of the quasi-singular 


.. /^2 y ,2 \ 

operator Vs = , since the hrst term in the righ-hand side of Eq. flTTl) is a regular 

expectation value which representes the Coulomb electron-electron repulsion. 

In this study we have calculated the expectation values of the Vs-operator for each of 
the two-electron ions considered in this study. These expectation values (in atomic units) 
can be found in Table III. It is interesting to investigate numerical changes of the () 

’’31 

/ 2 2 

and/or ^((^) — (^) expectation values in the series of two-electron ions, i.e. when the 
nuclear electric charge Q increases. By using our results from Table III one can solve this 
problem for the expectation values of the Vg operator and after a few simple recalculations 
for the expectation values of the Vitni operator. 


V. SCATTERING OF FAST ELECTRONS BY THE TWO-ELECTRON IONS IN 
THE BORN APPROXIMATION 

Another interesting problem which is of interest for the two-electron ions is the scat¬ 
tering of fast electrons at such ions. This problem is actual in physics of the hot atomic 
plasmas, which can be created in modern laboratory, used in various technical applications 
and observed in astrophysics. The problem is formulated as follows. As is well known the 
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electric charge density distribution p(r) in arbitrary atom/ion can be represented in the 
form p(r) = Qed(r) — en(r), where Qe is the electric charge of the central nucleus, ?7,(r) is 
the electron density and 5(r) is the electron-nucleus delta-function. The electron density is 
normalized by the condition / n(r)d^r = N^, where is the number of bounded electrons. 
For atoms we always have Nf> = Q, while in positively charged ions we have N^. < Q. For 
two-electron ions we have iVg = 2. From here one hnds the following expression for the Born 
amplitude of elastic scattering: 


/s(g) = fshk sin =-- J'e(g)] =-- f exp(-zq ■ r)n(r)d^r] (12) 

^ 2^ oog oog J 

where oq is the Bohr radius and the value is the electron form-factor of the ion (see, 


e.g., 


16|). For two-electron ions with large Q {Q > 2) the electron density is a spherically 


symmetric function, i.e. n{r) = n{r). In these cases we have 

= J exp(—zq ■ r)?7,(r)(i^r = 47r J ^ n{r)r‘^dr 


(13) 


For small values of the transfered momentum the kernel of the last integral in Eq. lIT^ is 
represented as a power series in qr. This leads to the following formula for 

2 


fB{q) = fB (2k sin 


aoq^ 


(Q _ AT,) + ^ - ^ + ^ - + 


3! 


5! 


7! 


9! 


(14) 


For very compact two-electron ions with Q >2 the series in the right-hand side of Eq. fll4p 
converges very fast. The differential cross-section in the Born approximation can easily be 
obtained from the formula, Eq. fll4p . Finally, we arive to the formula 


da 

dVt 


a^q 


q^{r'^) 


q,4(^4) 


+ ... 


(15) 


3!(Q-iV,) 5!(Q-iVe) 

where Q > N^. The hrst term in this formula corresponds to the usual Rutherford scattering, 
while other terms in Eq. flTsp represent the consequtive corrections for the electron scattering 
forward which is markedly anysotropic, since q^ = 4fc^sin^|. Largest deviations from the 
Rutherford scattering can be found for the small angles (i.e. for 6^ ~ 1 and qRa ^ 1, where 
Ra is the radius of the atom), where the cross-section of the scattering forward can exceed 
the cross-section of the Rutherford scattering. Note that application of these formulas to 
the two-electron H“ ion is not apppropriate (it leads to an answer which is inaccurate and 
has no direct physical sense). 
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VI. CALCULATIONS AND FORMULAS FOR THE Q ^ EXPANSIONS 


In this Section we discuss a few important details of actual calculations of the bound 
state properties of two-electron ions considered above. First of all, let us note that the total 
number of bound state properties evaluated in this study to very high accuracy signihcantly 
larger (for each two-electron ion) than a few properties shown in Tables I - III. In particular, 
we determine the both electron-nucleus and electron-electron cusp values, i.e. the following 
ratios of the expectation values 






(<^(rejv)) 


in the case of the electron-nucleus cusp, and 


= 


(^(’^ee)aree) 


(16) 


(17) 


(5(ree)) 

for the electron-electron cusp. These two expectation values must coincide with the known 
values of these cusps, i.e., with the following numerical values (in atomic units) 

2 rUeMisi _ 1 / _ _ 1 \ -1 


= -Qe 


= -Q 




-q(i + m^')’ 


= 0.5 


( 18 ) 


rUe + Mj^ 

where Mj^ = is the nuclear mass which can be hnite (real), or inhnite for model atomic 
systems. 

The coincidence of these two expectation values, Eqs. (lT6|) - (ITT)) , with the predicted values, 
Eqs. fllSp . is a very effective test for the variational wave functions in any Coulomb few-body 
system. In reality, this test is often used to show the correctness of our expectation values of 
the two-particle delta-functions determined with the use of variatonal expansion, Eq.(|3]). For 
instance, for the AT^+, Ca^®+ and Co^®+ ions the computed electron-nucleus cusp values are: 
UeN = -13.000000009571, UeN = -20.000000015543 and UeN = -27.000000010287, respectively. 
Analogous electron-electron cusp values are = 0.50000000545, Uee = 0.50000000652 and 
z/ee = 0.50000000875, respectively. This indicates clearly that all expectation values of the 
two-particle delta-functions used in this study are highly accurate. 

Second, in our calculations of the lowest-order QED corrections we have used the Q~^ 
expansion to evaluate the numerical values of the reduced Bethe logarithm, i.e. for the 
value: 


Infco = 



N 

^2 ^3 ^ 

'Q 


(19) 


10 









In actual computations this series converges very fast, since the coefficients Ci,C 2 , ■ ■ ■ are 
signihcantly smaller than the hrst coefficient Co (this is the main advantage to apply the 
reduced Bethe logarithm Infco in actual calculations). Therefore, we can restrict ourselves 
by a relatively small N in Eq. flT^ . In this study we chose N = 6. By using the known values 
of the Bethe logarithm for the two-electron He atom (Ini^o = 4.2701621981) and H~ (Inii'o 
= 2.993044153), Li’*', Be^C B^+, ions (the same values of Init'o which are mentioned 
in the hrst reference from [l^) we have found the following values of the Ck coefficients in 
Eq.(Il9D: 


Co = 2.9841349233 , Ci = -0.0124321314 Ca = 0.0235421443 (20) 

^3 = -0.0008950491 , C 4 = 0.0068528122 , C 5 = -0.0081585463 

Now, by using these six cefficients we hnd for the ions In/co = 2.9836752649 ~ 

2.9836753 and Ini^o = 9.4214269146 ~ 9.4214269. Numerical evaluation of Bethe loga¬ 
rithms for other two-electron ions is performed analogously. Note that the formula, Eq. flT^ . 
can be applied to the ground l^S*-states in the two-electron ions/atoms only. For the excited 
states the series, Eq. dT^ . must contain an additional term. 

Third, numerical calculation of the expectation value of the Vitni operator is 

reduced to computations of the ‘singular’ part of this operator Vs = n(-F-(the 

^ V'’31 ’’31 / 

corresponding numerical values can be found in Table III) and electron-electron repulsion 
energy in atomic units (^) = (’^ 21 ^)- These expectation values can be found in Table IV 
which also contains the expectation values of single-electron kinetic energies (—|Vg) and 
triple delta-functions (^ 321 ) (all values expressed in atomic units). 

In order to represent the expectation values determined in this study (see Tables I - IV) 
we can derive relatively simple and convenient interpolation formulas (or asymptotic for¬ 
mulas) which represent the Q—deppendence (or deppendence) of different expectation 
values. Such interpolation formulas are of great interest in actual applications to few-body 
systems, since in reality it is difficult to use (and even keep) huge tables of highly accurate 
data. However, if we know the explicit formula for the Q—expansions of these expectation 
values, then we can operate with a restricted number of the expansion coefficients deter¬ 
mined separately for each expectation value. Let us consider the (^^ 2 /) = (ree) expectation 
values determined for different two-electron ions. Numerical values of this expectation value 
for different ions can be found in Table IV. As follows from this Table the expectation value 
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(rgg^) increases with Q. Therefore, we can approximate these expectation values by using 
the following Q“^-expansion: 


{fee) — A2Q^ + AiQ + ^0 + ‘^ + ^ + 7tI + -- - + 7^ + ■ ■ ■ 




Bn 


( 21 ) 


where the hrst three terms represent the regular part of the Laurent expansion (or series), 
while all terms with the negative powers of Q represent the principal part of the Laurent 
series. By using our data for the {r~^) expectation values from Table IV we can estimate 
the hrst few coefficients ^ 2 , Ai, Aq, i?i, i? 2 ,-S 3 , • • • from a least-square ht to our calculated 
expectation values of up to Q = 28. The four different sets of such coefficients are 

shown in Table V. They correspond to the cases when the expansion Eq.(j2T]) contains 10, 12 
and 14 terms. It is intereting to illustrate the power of the general theory which is behind 


the Q“^-expansions (see, e.g., 15|]). As folows from the comparison of the hrst two columns 
in Table V our guess about the explicit form of Eq. (l21l) is wrong. Indeedn the hrst coefficient 
A 2 in the hrst column is almost zero, while the secod coefficient is signihcantly larger (in 10® 
times!). This means we have to assume that y42 = 0 identically, or, in other words, the hrst 
actual term in Eq. (l2T]) must be AiQ. This was corrected in calculations performed for the 
columns 2 - 4 in Table V. Briehy, this means that one can make a mistake in Eq. fl2T]) . but 
general theory will correct it. Briehy, we can say that the correct form of the q~^ expansion 
for each expectation value will be restored in any case, if the original computational data 
are truly highly accurate. 

Interpolation (or asymptotic) formulas for other expectation values can be derived analo¬ 
gously. Note that in actual applications one hnds three following situations: (1) expectation 
value decreases when the nuclear charge Q grows, (2) expectation value increases when the 
nuclear charge Q increases, and (3) expectation value is alomst constant during Q variations. 
For instance, all (rg^) and (r^^) expectation values with positive k represent the hrst group. 
The expectation values (rg^y) and (rgg) with negative k belong to the second group. The 
reduced Bethe logarithm mentioned above corresponds to the third group. The explicit for¬ 
mulas of the Q—expansions in each of these cases can be derived analogously. Furthermore, 
this procedure can be applied to each of the expectation values presented in Tables I - IV. 
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VII. CONCLUSIONS 


We have investigated the Q—dependencies of different bound state properties of the 
ground 1 ^5'—states in the seires of two-electron ions which includes all ions from the hydrogen 
two-electron ion up to the two-electron nickel ion. A substantial part of this series (ions 
from N®’*' to is of great interest, since these ions were not considered in earlier highly 

accurate studies. In particular, the lowest-order QED corrections for each of these ions 
have been determined to high accuracy for the first time. For convenience some results of 
such calculations are shown in Table VI where all results are given in atomic units. We also 
determine the specihc component of the isotopic shifts in these ions by using the expectation 
values of the quasi-singular (or Frullanian) Vinti operator. Scattering of fast electrons by the 
two-electron ions is briefly discussed in the Born approximation. A brief theory of accurate 
interpolation formulas for Q—dependencies of different bound state properties is discussed 
and illustrated with a few examples. 
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TABLE I: The total non-relativistic energies E and {r~^) and {r~g) expectation value for a number 
of two-electron atoms/ions in their ground 1^S'—states (in atomic units). All nuclear masses are 
assumed to be infinite. 


ion/atom 

Q 

E 

{’’' en ) 


H- 

1 

- 0.5277510165443771965925 

1.1166628245254385975 

0.1551041525624275657 

He 

2 

- 2.90372437703411959831115924519440 

6.0174088670242866935 

1.4647709233190671890 

Li+ 

3 

- 7.27991341266930596491875 

14.927623721321274552 

4.0822327875220922168 

Be2+ 

4 

- 13.65556623842358670208051 

27.840105671230351970 

8.0288017817568097761 

B3+ 

5 

- 22.03097158024278154165469 

44.753485397059351468 

13.307121456650923770 

C4+ 

6 

- 32.40624660189853031055685 

65.667312718097517695 

19.918006872526419836 

N5+ 

7 

- 44.781445148772704645183 

90.58139550837612183 

27.86179483571140529 

06 + 

8 

- 59.156595122757925558542 

119.4956378746231599 

37.13865011529274487 

F7+ 

9 

- 75.531712363959491104856 

152.4099865906842945 

47.74866291286534100 

Ne®+ 

10 

- 93.906806515037549421417 

189.3244097370080602 

59.69188677567529073 

Na®+ 

11 

- 114.28188377607272189582 

230.2388870075436700 

72.96835550416178458 

Mg“+ 

12 

- 136.65694831264692990427 

275.1534048677168016 

87.57809148852457738 


13 

- 161.03200302605835987252 

324.0679539487669773 

103.5211101497819791 

Sil2+ 

14 

- 187.40704999866292631487 

376.9825275594010858 

120.7974224560008489 

pl3+ 

15 

- 215.78209076353716023462 

433.8971207930373478 

139.4070364215928324 

S14+ 

16 

- 246.15712647425473932009 

494.8117299699472423 

159.3499580389784688 

C115+ 

17 

- 278.53215801540009570337 

559.7263522763435335 

180.6261918784910823 

Ari6+ 

18 

- 312.90718607661114879880 

628.6409855238003586 

203.2357414866408062 

K17+ 

19 

- 349.28221120345316700447 

701.5556279846605630 

227.1786096576231830 

Cai®+ 

20 

- 387.65723383315855621790 

778.4702782768293685 

252.4547986227706891 

Sci9+ 

21 

- 428.03225432023469116264 

859.3849352814886391 

279.0643101854956014 

Xi20+ 

22 

- 470.40727295513838395930 

944.2995980832556961 

307.0071458191845422 

V21+ 

23 

- 514.78228997811177388135 

1033.214265925954510 

336.2833067393932757 

01-22+ 

24 

- 561.15730558958127234352 

1126.128938179445216 

366.8927939578833437 

Mn23+ 

25 

- 609.53231995807574620568 

1223.043614314415181 

398.8356083236148089 

Fe24+ 

26 

- 659.90733322632780520901 

1323.958293882988089 

432.1117505542265143 

Co25+ 

27 

- 712.28234551602655145614 

1428.872976503642391 

466.7212212604835555 

Ni 26 + 

28 

- 766.65735693155709991040 

1537.787661849361804 

502.6640209654591693 
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TABLE II: The ((5(reAr)), (- 3 —)_r, ((5(ree)) and expectation values for a number of two-electron 

'’’eN ’’ee 

atoms/ions in their ground 1^5—states (in atomic units). The subscript R is used below to designate 
the regular part of the expectation value. All nuclear masses are assumed to be infinite. 


ion/atom 

Q 

{S{reN)) 

{^)r 

('5(ree)) 


H- 

1 

0.16455287284772 

-3.435594850575 

0.0027379921262 

0.064307887281 

He 

2 

1.81042931850139 

-53.67642660259 

0.106345370634 

-0.347101795446 

Li+ 

3 

6.85200943432227 

-238.7091644015 

0.533722536561 

-6.528108292960 

Be2+ 

4 

17.1981725446016 

-662.2679562154 

1.522895351487 

-26.72596510717 

B3+ 

5 

34.7587436609116 

-1437.181953992 

3.312442112829 

-70.59566341481 

C4+ 

6 

61.4435780565873 

-2682.795781398 

6.141043971072 

-148.7264629873 

N5+ 

7 

99.1625346267485 

-4523.555532659 

10.247411460455 

-272.4138789959 

06+ 

8 

149.825472814519 

-7088.104852970 

15.870266106959 

-453.5227647586 

F7+ 

9 

215.342252155782 

-10508.65528876 

23.248334282439 

-704.3954157123 

Ne®+ 

10 

297.622732139792 

-14920.52185328 

32.620344790203 

-1037.785453519 

Na9+ 

11 

398.576772218307 

-20461.76606779 

44.225027751185 

-1466.807903211 

Mgi°+ 

12 

520.114231810706 

-27272.91291287 

58.301114053422 

-2004.900136214 

Al“+ 

13 

664.144970311620 

-35496.72082712 

75.087335056216 

-2665.790465008 

Sii2+ 

14 

832.578847097505 

-45277.99108426 

94.822422421253 

-3463.472333880 

pl3+ 

15 

1027.32572153134 

-56763.40721339 

117.74510801161 

-4412.182728211 

S14+ 

16 

1250.29545296610 

-70101.39786688 

144.09412382866 

-5526.383843605 

CII5+ 

17 

1503.39790074724 

-85442.01834226 

174.10820197134 

-6820.747326894 

Ai.16+ 

18 

1788.54292421458 

-102936.8471892 

208.02607460887 

-8310.140582248 

K17+ 

19 

2107.64038270360 

-122738.8951897 

246.08647396232 

-10009.61476079 

Cai«+ 

20 

2462.60013554641 

-145002.5246139 

288.52813229153 

-11934.39414078 

Scl9+ 

21 

2855.33204207246 

-169883.3771019 

335.58978188593 

-14099.86666982 

Ti20+ 

22 

3287.74596160910 

-197538.3088595 

387.51015505776 

-16521.57548801 

V21+ 

23 

3761.75175348198 

-228125.3321072 

444.52798413712 

-19215.21128676 

Cr22+ 

24 

4279.25927701537 

-261803.5619184 

506.88200146822 

-22196.60538551 

Mn2^+ 

25 

4842.17839153237 

-298733.1677353 

574.81093940663 

-25481.72342952 

Fe24+ 

26 

5452.41895635514 

-339075.3289712 

648.55353031710 

-29086.65962870 

Co25+ 

27 

6111.89083080505 

-382992.1942021 

728.34850657198 

-33027.63147071 

Ni26+ 

28 

6822.50387420288 

-430646.8435300 

814.43460055003 

-37320.97485193 
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TABLE III: The {tIjj), {r‘tj<j) and Vs = \{(^) — (^) expectation values of the different two- 

V ’’31 ^31 

electron atoms/ions in their ground S'—states (in atomic units). All nuclear masses are assumed 
to be infinite. 


ion/atom 

Q 


(aV) 


H- 

1 

11.913699678051262 

645.14454241221937 

-0.4642618530805959 

He 

2 

1.1934829950189353 

3.9735649316629101 

-1.0010782750156071 

Li+ 

3 

0.44627901120132874 

0.52960197172324301 

-1.5089235297645549 

Be2+ 

4 

0.23206731553080959 

0.14056152936071842 

-2.0126039740505242 

B3+ 

5 

0.14196918178327382 

5.2083409587888660-10-2 

-2.5147510387558219 

C4+ 

6 

9.5739522935058123-10-2 

2.3539482331106725-10-2 

-3.0161594174263226 

N5+ 

7 

6.8904525396724106-10-2 

1.2141346213809833-10-2 

-3.5171547716283449 

06+ 

8 

5.1955356759041028-10-2 

6.8816453528063574-10-3 

-4.0178956847309996 

F7+ 

9 

4.0570590003835947-10-2 

4.1863570102541125-10-3 

-4.5184687145246674 

Ne®+ 

10 

3.2556160988739701-10-2 

2.6907921048456992-10-3 

-5.0189251377057326 

Na®+ 

11 

2.6701867132630897-10-2 

1.8073903567941317-10-3 

-5.5192972710023712 

Mg“+ 

12 

2.2295730542814468-10-2 

1.2585800135898979-10-3 

-6.0196064958067383 

Al“+ 

13 

1.889667256146306-10-2 

9.031562973287720-10-3 

-6.5198675240121025 

Sii2+ 

14 

1.621960949319265-10-2 

6.648054655635712-10-3 

-7.0200908109741794 

pl3+ 

15 

1.407363963658752-10-2 

5.001508016520978-10-3 

-7.5202839913624110 

S14+ 

16 

1.232702890357691-10-2 

3.834614539301480-10-3 

-8.0204527704226906 

CII 5 + 

17 

1.088649409859439-10-2 

2.989044580890498-10-3 

-8.5206014973391348 

Ari6+ 

18 

9.684480583887134-10-3 

2.364226371363055-10-3 

-9.0207335455792636 

K17+ 

19 

8.671093473064351-10-3 

1.894471957484780-10-3 

-9.5208515720108625 

Cai«+ 

20 

7.808833943501406-10-3 

1.535808715634216-10-3 

-10.020957697600605 

Scl9+ 

21 

7.069079045780357-10-3 

1.258148786675343-10-3 

-10.521053636055015 

Ti20+ 

22 

6.429674915299672-10-3 

1.040496858269455-10-3 

-11.021140787100487 

V21 + 

23 

5.873258366855471-10-3 

8.679410260860985-10-3 

-11.521220305246443 

Cr22+ 

24 

5.386065414197868-10-3 

7.297190978578113-10-3 

-12.021293151234381 

Mn2^+ 

25 

4.957071712852294-10-3 

6.179497295229088-10-3 

-12.521360131054423 

Fe24+ 

26 

4.577363304001367-10-3 

5.267838395821848-10-3 

-13.021421925898852 

Co25+ 

27 

4.239669745906067-10-3 

4.518271622832327-10-3 

-13.521479115417290 

Ni26+ 

28 

3.938013446919951-10-3 

3.897410055792940-10-3 

-14.021532195958377 
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TABLE IV: The (—^Vg) and ((532i) expectation values determined for the different two- 

electron atoms/ions in their ground states (in atomic units). All nuclear masses are assumed 
to be infinite. 


ion/atom 

Q 

(Ae') 


(< 5321 ) 

H- 

1 

0.3110215022143000515 

0.26387550827218860 

0.00512686786 

He 

2 

0.9458184487999231936 

1.45186218851705980 

1.86179097423 

Li+ 

3 

1.5677195591374732721 

3.63995670633465298 

33.3211589496 

Be2+ 

4 

2.1908707739064748156 

6.82778311921179335 

231.038802708 

B3+ 

5 

2.8146960470665581690 

11.0154857901213908 

996.009538154 

C4+ 

6 

3.4388907009922270508 

16.2031233009492652 

3221.84957561 

N5+ 

7 

4.0633059371614058208 

22.3907225743863523 

8596.07234611 

06 + 

8 

4.6878626272933908613 

29.5782975613789628 

19974.3019384 

F7+ 

9 

5.3125152146355572467 

37.7658561819797455 

41827.4545255 

Ne*+ 

10 

5.9372357244481370617 

46.9534032575187747 

80761.8501364 

Na®+ 

11 

6.5620060637014121725 

57.1409418880363610 

14611.2288400 

Mg“+ 

12 

7.1868140253485781079 

68.3284741563234650 

250608.071842 

Al“+ 

13 

7.8116510805680880023 

80.5160015130291799 

411111.980219 

Sil2+ 

14 

8.4365110930811231600 

93.7035249993314632 

49432.1960573 

pl3+ 

15 

9.0613895365056378705 

107.891045381768580 

993207.180907 

S14+ 

16 

9.6862829995969480726 

123.078563237127370 

1476863.50256 

C115+ 

17 

10.311188863149912371 

139.2660790077000478 

2142646.61319 

Ai.16+ 

18 

10.936105082961410952 

156.4535930383055744 

3041724.57849 

K17+ 

19 

11.561030040385195148 

174.6411056017265835 

4235364.75740 

Cai®+ 

20 

12.185962437154670069 

193.8286169165792781 

5796183.43295 

Scl9+ 

21 

12.810901219907308976 

214.0161271601173455 

7809468.39351 

Ti20+ 

22 

13.435845525071932709 

235.2036364775691919 

10374574.4651 

V21 + 

23 

14.060794637989284133 

257.3911449890558868 

13606391.9942 

CI. 22 + 

24 

14.685747962156786948 

280.5786527947906359 

17636888.2813 

Mn2®+ 

25 

15.310704995789691498 

304.7661599790378727 

22616721.9649 

Fe24+ 

26 

15.935665313746399507 

329.9536666131639021 

28716930.3559 

Co25+ 

27 

16.560628553438948983 

356.1411727580132750 

36130689.7235 

Ni26+ 

28 

17.185594403740284234 

383.3286784657785490 

45075148.5302 
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TABLE V: The coefficients of the Q~^ expansion determined for the (r“^) expectation values 
determiend for all 28 ions presented in Table IV (in atomic units). N is the total number of terms 
used in the Q~^ expansion. 


cefficient 

N = 10 

N = 10 

N = 12 

W = 14 

^2 

-0.494831796599196E-08 

O.OOOOOOOOOOOOOOOE+00 

O.OOOOOOOOOOOOOOOE+00 

O.OOOOOOOOOOOOOOOE+00 


0.625000467182303E+00 

0.624999946424256E+00 

0.625001683616873E+00 

0.625003246752918E+00 

^0 

-0.315350485009514E+00 

-0.315327955933086E+00 

-0.315514341765753E+00 

-0.315701190374814E+00 

Bi 

0.264528972963422E-01 

0.259297624553894E-01 

0.343629291054843E-01 

0.439685640740424E-01 

Bi 

-0.764173015283037E-02 

-0.468906451156542E-03 

-0.212091717363386E+00 

-0.492221452096754E+00 

Bi 

0.217453679005328E-01 

-0.386848769404071E-01 

0.322514909365720E+01 

0.837712175158511E+01 

Bi 

-0.104596000674926E+00 

0.211498434740675E+00 

-0.322114782408164E+02 

-0.950787058977685E+02 

Bs 

0.203688490986952E+00 

-0.807975095629530E+00 

0.210536612541492E+03 

0.731912378303186E+03 

Be 

-0.206581093655412E+00 

0.168265813518147E+01 

-0.900343198369528E+03 

-0.385843989156749E+04 

Br 

0.683035978612077E-01 

-0.177678156292532E+01 

0.246014385082417E+04 

0.138443324014201E+05 

Bs 


0.705173625864056E+00 

-0.406809330705356E+04 

-0.330327568745938E+05 

B9 



0.362963884774198E+04 

0.499600243794525E+05 

Bio 



-0.130271721358420E+04 

-0.437013456795518E+05 

Bn 




0.182548118898741E+05 

Bi2 




-0.211138704685216E+04 
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TABLE VI: The lowest order QED-corrections Eq-Q, and numerical values of 

In Aq, in atomic units. In this Table the notations AE^^^ stands for the sum of 

^ee ^ee ’ 

the last two terms from Eq.([6]) multiplied by a factor ^ = M (nuclear mass). All nuclear masses 
are assumed to be infinite. 


ion/atom 

Q 

InKo 

' pp 

' pp 

^-^OO 


H- 

1 

2.9930441530 

0.987145110780E-01 

-0.1367762464713E+01 

0.12485643657E-05 

0.555264710871E-07 

He 

2 

4.37016021981 

0.989273545062E+00 

-0.3092590081520E+02 

0.22261831904E-04 

-0.135298475377E-04 

Li+ 

3 

5.179849129 

0.178846906701E+00 

-0.1526042743967E+03 

0.11024751764E-03 

-0.197385810727E-03 

Be^+ 

4 

5.755091813 

-0.758869771350E+01 

-0.4461493461302E+03 

0.33035865562E-03 

-0.116481400438E-02 

B3+ 

5 

6.201467201 

-0.289702881864E+02 

-0.1000390699060E+04 

0.75812296982E-03 

-0.440940827548E-02 

C4+ 

6 

6.566235883 

-0.715558284877E+02 

-0.1910673007667E+04 

0.14744207461E-02 

-0.127893001871E-01 

N5+ 

7 

6.874659414 

-0.143641108746E+03 

-0.3277442371480E+04 

0.25617115664E-02 

-0.310542023859E-01 

06 + 

8 

7.141831512 

-0.254091119110E+03 

-0.5205342434111E+04 

0.41010004202E-02 

-0.664058639717E-01 

F7+ 

9 

7.377493052 

-0.412248230953E+03 

-0.7802584739240E+04 

0.61692976733E-02 

-0.129086002178E+00 

Ne®+ 

10 

7.588297026 

-0.627866111317E+03 

-0.1118048429795E+05 

0.88374333231E-02 

-0.232987591610E+00 

Na®+ 

11 

7.778989583 

-0.911059814059E+03 

-0.1545310262982E+05 

0.12168138196E-01 

-0.396286477650E+00 

Mgio+ 

12 

7.953075479 

-0.127226672979E+04 

-0.2073696471413E+05 

0.16214333661E-01 

-0.642090983605E+00 

Ahi+ 

13 

8.1132164336 

-0.172221518425E+04 

-0.3466213164326E+06 

0.21017588805E-01 

-0.109131248829E+01 

Sii2+ 

14 

8.2614815231 

-0.227189863118E+04 

-0.3481549672596E+05 

0.26606715083E-01 

-0.150232170036E+01 

pl3+ 

15 

8.3995110079 

-0.293255406291E+04 

-0.4385365145496E+05 

0.32996475895E-01 

-0.219369069970E+01 

Sl4+ 

16 

8.5286271993 

-0.371564368022E+04 

-0.5438972182745E+05 

0.40186393648E-01 

-0.312285083765E+01 

C115+ 

17 

8.6499116644 

-0.463283913392E+04 

-0.6654976314062E+05 

0.48159640601E-01 

-0.434783445856E+01 

Ari6+ 

18 

8.7642603335 

-0.569600783126E+04 

-0.8046135394383E+05 

0.56882002451E-01 

-0.593579835268E+01 

K17+ 

19 

8.8724236914 

-0.691720092580E+04 

-0.9625350501886E+05 

0.66300905715E-01 

-0.796376203590E+01 

Cai«+ 

20 

8.9750366476 

-0.830864269774E+04 

-0.1140565786357E+06 

0.76344501477E-01 

-0.105193554360E+02 

Sci9+ 

21 

9.0726411126 

-0.988272109625E+04 

-0.1340022164342E+06 

0.86920799336E-01 

-0.137015754571E+02 

p^20+ 

22 

9.1657033175 

-0.116519792627E+05 

-0.1562232746201E+06 

0.97916846311E-01 

-0.176215509500E+02 

V21+ 

23 

9.2546272805 

-0.136291078896E+05 

-0.1808537654137E+06 

0.10919794629E+00 

-0.224033156641E+02 

0^22+ 

24 

9.3397654052 

-0.158269382973E+05 

-0.2080288038885E+06 

0.12060691617E+00 

-0.281845887975E+02 

Mn2^+ 

25 

9.4214269146 

-0.182584361318E+05 

-0.2378845594865E+06 

0.13196337542E+00 

-0.351175627958E+02 

Fe24+ 

26 

9.4998846295 

-0.209366956035E+05 

-0.2705582116209E+06 

0.14306306617E+00 

-0.433696980840E+02 

Co25+ 

27 

9.5753804665 

-0.238749342007E+05 

-0.3061879088677E+06 

0.15367720134E+00 

-0.531245244377E+02 

Ni26+ 

28 

9.6481299371 

-0.270864878203E+05 

-0.3449127313289E+06 

0.16355183851E+00 

-0.645824487293E+02 
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